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Abstract 

We consider Noether symmetry approach to find out exact cosmological solutions in f(T)- 
gravity. Instead of taking into account phenomenological models, we apply the Noether symmetry 
to the f(T) gravity. As a result, the presence of such symmetries selects viable models and allow 
to solve the equations of motion. We show that the generated f(T) leads to a power law expansion 



Of 



for the cosmological scale factor. 



1 Introduction 

The expansion of the universe is currently undergoing a period of acceleration which is directly mea- 
sured from observations such as Type la Supernovae [T], [2] cosmic microwave background (CMB) 
radiation [3j 0] , large scale structure [5] , baryon acoustic oscillations [6] , and weak lensing [7] . There 
are two remarkable approaches to explain the late time acceleration of the universe: One is to in- 
troduce some unknown matters called "dark energy" in the framework of general relativity (for a 
review on dark energy, see, e.g., [9]). The other is to modify the gravitational theory, e.g., f(R) 
gravity [lOlfTT]. 

To extend gravity beyond general relativity, "teleparallelism" could be considered by using the 
Weitzenbock connection, which has no curvature but torsion, rather than the curvature defined by 
the Levi-Civita connection |12j . |13j . This approach was also taken by Einstein |14j . The teleparallel 
Lagrangian density described by the torsion scalar T has been promoted to a function of T, i.e., f(T), 
in order to account for the late time cosmic acceleration |15j as well as inflation |16| . This concept 
is similar to the idea of f{R) gravity. Various aspects of f(T) gravity have been examined in the 
literature |17[ [18} [l~9j [201 El]- In particular, the presence of extra degrees of freedom and the violation 
of local Lorentz invariance as well as the existence of non-trivial frames for f(T) gravity have been 
noted [21]. Evidently, more studies on f(T) gravity are needed to see if the theory is a viable one. 
For a comprehensive review of the teleparallel gravity, the reader is referred to |22| . 

In this paper, we consider a flat FRW space-time in the framework of the metric formalism of f(T) 
gravity. Following |23j, we set up an effective Lagrangian in which the scale factor a and torsion scalar 
T play the role of independent dynamical variables. This Lagrangian is constructed in such a way that 
its variation with respect to a and T yields the correct equations of motion as that of an action with 
a generic f(T) mentioned above. The form of the function f(T) appearing in the modified action is 
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then found by demanding that the Lagrangian admits the desired Noether symmetry |24[ [25] . For a 
study of the Noether symmetry in f(R) cosmology see |26j. Similarly, we shall see that by demanding 
the Noether symmetry as a feature of the Lagrangian of the cosmological model under consideration, 
we can obtain the explicit form of the function f(T). Since the existence of a symmetry results in 
a constant of motion, we can integrate the field equations which would then lead to a power law 
expansion of the universe. 



= VABe^e^ , (1) 



2 f{T) gravity and cosmology 

To consider teleparallelism, one employs the orthonormal tetrad components eA(x p ), where an index 
A runs over 0,1,2,3 to the tangent space at each point x^ 1 of the manifold. Their relation to the 
metric g^ v is given by 

where [i and v are coordinate indices on the manifold and also run over 0,1,2,3, and forms the 
tangent vector on the tangent space over which the metric t]ab is defined. 

Instead of using the torsionless Levi-Civita connection in General Relativity, we use the curva- 
tureless Weitzenbock connection in Teleparallelism |12] . whose non-null torsion T p fll/ and contorsion 
are defined by 

rV = <( 9 /-^) , (2) 

K^ p = -\(T^ p -T»%-T p ^) , (3) 

respectively. Moreover, instead of the Ricci scalar R for the Lagrangian density in general relativity, 
the teleparallel Lagrangian density is described by the torsion scalar T as follows 

r^ P %, (4) 

where 

Sp ixv = 1 ^ + 5 , ^ _ 5 u T a^ (fi) 

The modified teleparallel action for f(T) gravity is given by 

d 4 z|e|/(T), (6) 

where |e| = det (e4) = \/—g and the units has been chosen so that c = IQttG = 1. Note that in 
action ([6]), we have overlooked any matter contribution in the action. Varying the action in Eq. © 
with respect to the vierbein vector field ei, we obtain the equation |15| 

\d, {eS/ v ) h - e\T\ x S p + S/»d, (T) f TT + \e A f = 0, (7) 

where a subscript T denotes differentiation with respect to T. We assume the four-dimensional flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space-time with the metric, 

ds 2 = h a pdx a dxP + r 2 dn 2 , (8) 

where f = a(t)r, x° = t and x 1 = r with the two-dimensional metric h a p = diag( — 1, a 2 (t)). Here, a(t) 
is the scale factor and dti 2 is the metric of two-dimensional sphere with unit radius. In this space-time, 
g^y = diag( — 1, a 2 , a 2 , a 2 ) and the tetrad components = (1, a, a, a) yield the exact value of torsion 
scalar 

T = -6H 2 , (9) 
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where H = a/a is the Hubble parameter and the dot denotes the time derivative of d/dt. 

In the flat FLRW background, it follows from Eq. ([7]) that the modified Friedmann equations are 
given by |15j 

12f T H 2 + f = 0, (10) 



H 



1 



4T/ TT +2/ 7 



-Tf T + 



f 



(11) 



It is known that f(T) gravity has first-order gravitational field equation in derivatives. Similar to 
f(R) gravity in general relativity where the gravitational field equation is fourth-order in derivatives, 
it is important to investigate the theoretical aspects in order to examine whether f(T) gravity can be 
a gravitational theory like f(R) gravity. 

In order to derive the cosmological equations in a FLRW metric, one can define a canonical 
Lagrangian C = C(a,a,T,T), where Q = {a,T} is the configuration space and TQ = {a,a,T,T} 
is the related tangent bundle on which C is defined. The variable a(t) and T(t) are the scale factor 
and the torsion scalar in the FLRW metric, respectively. One can use the method of the Lagrange 
multipliers to set T as a constraint of the dynamics. Selecting the suitable Lagrange multiplier and 
integrating by parts, the Lagrangian C becomes canonical. In our case, we have 



2tt z 



dta 3 \f(T)-\ 



T + 6 



(12) 



where a is the scale factor scaled with respect to today's value (so that a = a/do an d ci(io) = !)• 
This choice for a, makes it dimensionless, whereas [/] = M 4 . It is straightforward to show that, for 
f(T) = —T, one obtains the usual Friedmann equations. 

The variation with respect to T of the action gives A = fx- Therefore the previous action can be 
rewritten as 



I = 2ir 2 dta 3 \f-f T 



T + 6 ^ 



and then, integrating by parts, the point-like FLRW Lagrangian is obtained 

£ = a 3 (f-f T T)-6f T ad 2 , 



(13) 



(14) 



which is a canonical function of two coupled fields, T and a, both depending on time t. The momenta 
conjugate to variables a and T are 

dC 

Va = — - = -12 f xaa. (15) 



da 



dC 
df 



Pt = — = 0. 

The total energy Ec corresponding to the 0, 0- Einstein equation is 

Ec = -6f T ad 2 -a 3 (/-/ T T) = 0. 
The equations of motion for T and a are 



a 3 /: 



■/■■/' i I 6 fj 
1 a z 



0, 



- 6/ T H 2 - 12/r- = 3(/ - f T T) + 12f TT f H, 

respectively. Note that from the equations ([17p . (|18|) and (|19|) we can recover the equation (|1U|) . 
and (jlip . respectively. Here we have taken Jtt 7^ 0. 



(16) 
(17) 

(18) 
(19) 
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Considering T and a as the variables, we see that T coincides with the definition of the torsion 
scalar in the FLRW metric (excluding the case f TT = 0). Geometrically, this is the Euler constraint 
of the dynamics. Furthermore, as we will show below, constraints on the form of the function f(T) 
can be achieved by asking for the existence of Noether symmetries. On the other hand, the existence 
of the Noether symmetries guarantees the reduction of dynamics and the eventual solvability of the 
system. 



3 The Noether symmetry 

Solutions for the dynamics given by the point-like Lagrangian (I14p can be obtained by selecting 
cyclic variables related to some Noether symmetry |26[ I27| . In principle, as a physical criterion, this 
approach allows one to select f(T) gravity models which are compatible with the Noether symmetry. 
Let C(q l ,q l ) be a canonical, non degenerate point-like Lagrangian so that 



— = 0, detHij 



d 2 c 



dq L dqi 



¥= 0, (20) 



where i?y is the Hessian matrix related to the Lagrangian C and a dot denotes derivative with respect 
to the affine parameter U, namely the cosmic time t. The Lagrangian in analytic mechanics is of the 
form 

£ = T(q,q)-y(q), (21) 

where T and V are the positive definite quadratic 'kinetic energy' and 'potential energy', respectively. 
The energy function associated with C is: 

*"w- £ - (22) 

which is the total energy T + V and the constant of motion. Since our cosmological problem has 
a finite number of degrees of freedom, we are going to consider only point transformations. Any 
invertible transformation of the 'generalized positions' Q l = Q l (q) induces a a local transformation of 
the 'generalized velocities' such that: 

Q*(q) = f£<* (23) 

the matrix J = \\dQ i /dq^\\ is the Jacobian of the transformation on the positions, and it is assumed 
to be non-zero. A point transformation Q l = Q l (q) can depend on a (or more than one) parameter. 
As starting point, we can assume that a point transformation depends on a parameter e, so that 
Q 1 = Q l (l> e )> an d that it gives rise to a one-parameter Lie group. For infinitesimal values of e, the 
transformation is then generated by a vector field: for instance, d/dx is a translation along the x axis 
while x{d/dy) — y{d/dx) is a rotation around the z axis and so on. 

In general, an infinitesimal point transformation is represented by a generic vector field on Q 

x = « l (q)|i- ( 24 ) 



The induced transformation (|23j) is then represented by 



A function -F(q, q) is invariant under the transformation X if 
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where LxF is the Lie derivative of F. If LxC = 0, X is then a symmetry for the dynamics derived 
by C Let us now consider a Lagrangian C and its Euler-Lagrange equations 

dXdqi dqi ' V ' 

Contracting (|27p with a l s leads to 

d dC dC 



Using the total derivative relation as 



dXdqi dqi< °' (28) 



, d dC d ( idC\ fda>\ dC 



dX dqi dX \ dqi J \ dX J dqi ' 
we obtain from equation (|28|) that 

T>{ a 'w) =LxC - (30) 

The immediate consequence is the Noether theorem which states: if Lx£- = 0, then the function 

^o = « fc ^, (31) 

is a constant of motion. Equation (|3ip can be expressed independently of coordinates as a contraction 
of X by a Cartan 1-form 

0C - §dq>. (32) 

For a generic vector field Y = y l d/dx l , and 1-form j3 = fiidx 1 , we have by definition iy/3 = y i Pi. Thus 
equation (|3ip can be expressed as: 

ixOc = S . (33) 
Through a point transformation, the vector field X becomes: 

* = UxdQ k ) + ( 4^(ixdQ k )) -2-. (34) 



If X is a symmetry and we choose a point transformations such that 

ixdQ 1 = 1; = i^l, (35) 

we get 

*=a§r; W = °- < 36 » 

Consequently Q 1 is a cyclic coordinate, related to conserved quantities, reducing the dynamics of the 
system to a manageable one. Furthermore, the change of coordinates given by (|35p is not unique. In 
general, the solution of equation (I35p is not defined on the whole space, rather, it is local in the sense 
explained above. 

Following |24[ I25| , we define the Noether symmetry induced on the model by a vector field X on 
the tangent space TQ = ^a,T,d,T^j of the configuration space Q = (a,T) of Lagrangian (|14p 

d „ d da d dB d ,„„, 

such that the Lie derivative of the Lagrangian with respect to this vector field vanishes 

L X C = 0. (38) 
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In ([37]) . a and /3 are functions of a and T and -jt represents the Lie derivative along the dynamical 
vector field, that is, 

d d ■ d , . 

(it da dT y ' 

It is easy to find the constants of motion corresponding to such a symmetry. Indeed, equation (I38p 
can be rewritten as 

r „ f dC dadC\ ( BL d/3 dC\ n 
Noting that ^ = we have 

dpa .da \ { q dpT dp \ 
v «— + — Pa ) + ( /3— + — PT ) = 0, (41) 

which yields 



dt dt J V (it cit 



^ (ap a + /? PT ) = 0. (42) 



Thus, the constants of motion are 



Q = ap a + Ppr, (43) 
whereas in f(T) theory we have px = and so the corresponding constant of motion becomes Q = ap a . 



In order to obtain the functions a and j3 we use equation (l40j) . In general, this equation gives a 
polynomial in terms of a 2 and dT with coefficients being partial derivatives of a and (3 with respect to 
the configuration variables a and T. Thus, the resulting expression is identically equal to zero if and 
only if these coefficients are zero. This leads to a system of partial differential equations for a and f3. 



4 Noether symmetries in f(T) cosmology 

For the existence of a symmetry, we can write the following system of equations (linear in a and /3) , 

f T (a + 2ad a a) + af TT (3 = 0, (44) 

af T d T a = 0. (45) 

which are obtained by setting the coefficients of the terms d 2 and dT in Ly^C = to zero. In order 
to make L-^C = vanish we will also look for those particular f's which, given the Euler dynamics, 
also satisfy the constraint 

3a(f-Tf T )-a(3Tf TT = 0. (46) 

This procedure is different from the usual Noether symmetry approach, in the sense that now L-^C = 
will be solved not for all dynamics (which solve the Euler-Lagrange equations), but only for those / 
which allows Euler solutions to solve also the constraint (I46h . Imposing such a constraint on the form 
of / will turn out to be, as we will show, a sufficient condition to find solutions of the Euler-Lagrange 
equation which also possess a constant of motion, i.e. a Noether charge. As we shall see later, the 
system f|44[) and (j45]) can be solved exactly. Having a non-trivial solution for a and j3 for this system, 
one finds a constant of motion if also the constraint (I46p is satisfied. A solution of (|44D and (l45l) exists 
if explicit forms of a, /3 are found. If, at least one of them is different from zero, a Noether symmetry 
exists. 

When from equation (l45j) we get 

d T a = => a = a(a). (47) 

Thus, from this equation it can be seen a only depend on a. On the other hand, we can obtain a 
from equation (j46|) as follows 

«W = 3(7^T)«"' T) ' (48) 
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If one uses this expression in equation (|44p to eliminate a(a), one obtains 

To solve this equation we assume that the function f3(a, T) can be written in the form /3(a, T) = 
A(a)B(T), where A and are separate functions of a and T, respectively. Substituting this ansatz 
for j3(a,T) into equation f|49[> . we obtain 

2a cL4 3/ 



^4 da T/ T 



(50) 



Since the left-hand side of this equation is a function of a only while the right-hand side is a function 
of T, we should have 

-^L = C = Const, (51) 

which results in 

f(T) = f Q T(-v). (52) 
On the other hand equation (|50p . with its right hand-side equal to C, has the solution 

2 ^ = C ^A(a) = a c / 2 . (53) 

A da w v ' 

Now, using ([EI]) and /3(a,T) = a c / 2 B{T) in equation (JIB"]) , we find 

a(o) = ^T- 1 a§ +1 J B(T). (54) 
Since a(a) should be a function of a only, from the above expression for a(a) we can write 

B(T) = T, (55) 

and thus we obtain 

(3(a,T) = Ta c/2 , a(a) = ^a§ +1 . (56) 

To obtain the corresponding cosmology resulting from this type of f(T), we note that the existence 

of Noether symmetry implies the existence of a constant of motion Q = ap a . Hence, using equation 
(|15p we have 

q = M .'frw+»>. (57) 



On the other hand, by taking f TT ^ from equation (j 1 8[) we obtain 

T = -6^. (58) 

By using equations (|57p and (|58p we can obtain scale factor as following 

o(t)~ti?. (59) 

Therefore, in the context of f(T) = f T~^ cosmology, the universe evolves with a power law 
expansion. It is seen from (f59|) that the condition under which the universe expands is C < 0. The 
deceleration parameter as a function of C is therefore given by 

q(C) = - (i + y ) (60) 

The condition for acceleration is q(C) < 0, thus we have C > — 2. As it can be seen, for C — > we have 
q — > — 1, that is the universe finally approaches the eternal de Sitter phase with infinite acceleration, 
however, the accelerated expansion occurs in — 2 < C < 0. 
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Figure 1: The behavior of effective EOS, function of z. When C 

behavior in the model. 



-1 and f2 m = 0.33, we get a phantom 



5 Dark energy equation of state and age of the Universe in f(T) 

In this section, we consider the accelerated expansion of the universe in the context of f(T) theory, 
without introducing the mysteries fluid the so-called "dark energy" with a negative equation of state 
(EOS) parameter w. Let us then start from equation (1591) and find the Hubble parameter as a function 
of the redshift z as 

9 n 

(61) 



H(z) = --H (l + z) 



c 

2 



where ao/a = 1 + z with ao and Hq being the values of the parameter at the present epoch. It is worth 
noting that equation (16 1 P is the same as that derived from the standard Friedmann equation with 
w = —2/3. Now, we may write the Friedmann equation in a formal fashion which would encapsulate 
any modification to the standard Friedmann equation in the last term regardless of its nature |28| 
that is 



H 2 /Hl = n m (l + zf + 5H 2 /H { 



o ■ 



where VL m = p/po c , Po c = 3Hq. Also, defining the effective EOS, denoted by w cS (z), as 



1 dln5H 2 

+ 



3dln(l + z)' 

we can calculate w eB (z) using equations ([6ip . ([62p and (|63p with the result 



-1 + 



A 

(71 



[i + z y 



O m (l + z) 



(62) 



(63) 



(64) 



Figure 1 shows the behavior of the effective EOS parameter, function of z. As it can be 

■1, which is the characteristic of one 



0.33 we have w „ ^ 

en r~Z> 



seen, for C = — 1 (-2 < C < -1) and ft r , 
type of dark energy, the so-called phantom and from equation (|64p . for C — > 0~, we have w oS — > — 1. 

To continue we consider the age of the Universe in f(T) model. Thus, the age of the matter 
dominated Universe in FLRW models is given by 



Hoto 



1 



3 VI - 



In 



(65) 



where Hq 1 = 9.8 x lO 9 ^" 1 years and the dimensionless parameter h, according to resent data, is 
about 0.7. Hence, in the flat matter dominated universe with Vt tota i = 1 the age of the universe would 
be only 9.3 Gyr, whereas the oldest globular clusters yield an age of about ~ 13.5 Gyr [29j. We 
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Figure 2: Hoto as a function of Q m for C = —1.9 (solid line), C = —3 (dashed line) and C = —6 (dotted line). Figure 
shows that for a fixed value of fl m the predicted age of the universe is longer for larger values of C. 



obtain the age of the universe for our model from equations 
the Friedmann equations as follows 



TUJ), (|52P and pJJ by taking matter in 



Hoto 



1 



C ^T^W r 



In 



(66) 



For a flat, matter dominated Universe with Q m ~ 0.3 and C = —3 we have a prediction for the age of 
the Universe of about 13.2 Gyr. It seems that the age of the universe in our model is longer than the 
FLRW model. Figure 2 shows the behavior of the dimensionless age parameter, Hoto, as a function 
of O m for different values of C. 



Note added 

Before sending our present work to the Archive (on 13 December 2011), we became aware of the 
fact that there already existed a paper arXiv:1112.2270 [gr-qc] [30] (sent to the Archive on 10 De- 
cember 2011), in which the authors had addressed similar issues. Although the content and results 
of the present work overlap with the ones in |30j, it is worth noting that our work was completed 
independently. 
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6 Conclusions 

In this paper we have studied a generic f(T) cosmological model by using Noether symmetry approach. 
We have taken the background geometry as a flat FLRW metric and derived the general equations 
of motion in this background. The phase space has been constructed by taking the scale factor 
a and torsion scalar T as two independent dynamical variables. The Lagrangian of the model in 
the configuration space spanned by {a,T} is so constructed that its variation with respect to these 
dynamical variables yields the correct field equations. The existence of Noether symmetry implies 
that the Lie derivative of this Lagrangian with respect to the infinitesimal generator of the desired 
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symmetry vanishes. By applying this condition to the Lagrangian of the model, we have obtained the 
explicit form of the corresponding f(T) function. We have shown that this form of f(T) results in a 
power law expansion for the scale factor and the accelerated expansion occurs in —2 < C < 
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